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Abstract 

A new formulation for the study of interacting bosons on a lattice is in- 
troduced. This approach is used to give analytical expressions for the Mott 
insulating lobes in the phase diagram and to calculate the density-density 
correlation function. It is also shown that, at mean-field level, this newly 
introduced slave boson theory coincides with mean-field theory of a suitably 

introduced order parameter. 
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I. INTRODUCTION 



There has been recently a revival of interest for interaction-induced metal-to-insulator 
transitions. This is due to a large extent to the discovery of high-T c materials fl[] as in the 
cuprates it is believed that the insulating state originates from a strong Coulomb interaction 
in a half-filled band. Anderson's proposal that the Hubbard Model should capture the 
essential physics of the cuprates generated the development and the application of numerous 
theoretical techniques to the problem of interacting fermions on the lattice. In parallel these 
techniques have been applied to interacting bosons as well. In a similar way a Mott-insulating 
state should appear at commensurate densities due to a strong local interaction that plays 
the role of the Pauli principle and prevents the bosons from undergoing a Bose-Einstein 
condensation as known from the ideal Bose gas. This scenario has been recently discussed 
by Fisher et al || in the framework of a suitable mean-field theory. Their result for the 
pure case has been confirmed by Krauth et al. using Gutzwiller type of wave-function, 
and Sheshadri et al 0. Other mean-field theories have been proposed 0,0], and Quantum 
Monte-Carlo simulations all confirming Fisher et al.'s result. An alternative technique 
is provided by introducing auxiliary bosonic fields. The slave boson representations were 
pioneered in the context of spin models by Holstein and Primakoff and, in particular, by 
Schwinger who introduced a two slave boson representation of spin 1/2 ||10|| . A second class 
of slave boson representations was introduced by Barnes and employed by others, in the 
context of the Anderson model of a magnetic impurity in a metal |T| . The method has been 



extended to lattice fermion models fll2|-|15||, and gave rise to an extensive literature. In this 
letter I introduce a new representation for both interacting bosons and spins on the lattice 
and I apply it to the Bose-Hubbard model. I point out a close connection between slave 
boson mean-field theory and the collective field method. The density-density correlation 
function in the first Mott insulating lobe is calculated. 
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II. FORMULATION 



The new representation for both interacting bosons and spins on a lattice is based on a 
generalization of the slave boson representation of the fermionic Hubbard Model introduced 



by Kotliar and Ruckenstein JT2[ , where each atomic state is obtained with the help of distinct 
slave particles subject to constraints. Here the normalized n-fold occupied state at a given 
site i on the lattice is obtained as: 

| n >i= b^ A | vac > n > 0. (1) 

Even the empty site | >j is constructed by operating with the Bose creation operator 
onto a new vacuum state | vac >j meaning that even the empty lattice site is not pre-existent 
but rather created out of a total vacuum | >«= 6^ | vac >, in the enlarged Hilbert space. 
The slave bosons b n ^ are subject to the constraint: 

E KM = 1 ( 2 ) 

n>0 

indicating that, at each lattice site and at each time, the site has a well defined occupancy. 
The eigenvalues of the atomic problem with on-site interaction U 

E n = Un{n — 1) — fin (3) 

serve as chemical potential for the slave particles. In terms of those the physical boson 
creation operator reads: 

B+ = £ Vl+n- bt +hi b n ,i (4) 

n>0 

It is a straightforward exercise to show that B satisfies the canonical commutation relation 
provided the slave fields do it as well. This is sufficient in order to rewrite the original 
Hamiltonian. Taking as an example the Bose-Hubbard Model: 

h = y.h,n; />', + E(^ + ^ - 1) - v) B t^ (5) 

The corresponding partition sum as expressed in the slave boson formalism reads: 
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Z = J n^Db+^Dbn+^iDXi exp (-S) (6) 

with the action 

S = j P dTj2bt,i(dr + E n + i\i)b nti -J2 Ai 

" n,i i 

+ ^ Vl + n Vl + m b+ +lti b nti t itj b m+1 jb+ tj (7) 

where the A-field enforces the constraint. 

The limit of infinite local interaction simply results in restricting the number of slave 
fields to be 2. In this case it is worth noting that the b- fields can be taken as fermionic and 
the action (0) describes a quantum X — Y model in an external magnetic field p/pb- But 
in the following I shall keep the 6-fields as bosonic. The action ([/]) is invariant under a U(l) 
gauge transformation which allows to eliminate the phase of one slave field at the price of 
introducing a time-dependent A-field. 



III. RESULTS 

a) HARD-CORE LIMIT. 

As a first example of handling the action (|7|) I consider the hard-core limit in the saddle- 
point approximation. This yields, for | p, |<| to \, where to = t(k = 0): 

bl = *2±£ , b\ = ^ and A = ^ (8) 
2t ' 1 2t 2 K ' 

and both phase boundaries of the insulating state are properly recovered. Namely they read 
fx = —2dt and fi = 2dt for the empty (b\ = 0) and the full (bo = 0) systems. Whereas the 
density p is given by | b\ | 2 , the superfluid density as taken from 0] p s ~|< B >| 2 =| b bi | 2 is 
very different from p and they only coincide in the low density limit. To obtain the gaussian 
fluctuations I first take advantage of the U{1) gauge symmetry of S eq. (0) in order to take 
the &o-field as real. I can then integrate out both 6 and A fields. After having introduced 
the mean-field parameters the action reads: 
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S gf - r ^2(b[ b" _ k ) 
z fc 



Oik u n 

\ —u n (3k j 



(v \ 

°i,k 



(9) 



where oik = b 2 (tk — to) — Ab\tk and (3k = &o(£fc — *o)- This yields the spectrum 



2 " (t k - t ) 2 - 46^2 tfc ( tfc _ to ) = k ^ 2t 2 b 2 blk 2 (10) 



which is linear in k for small k as expected in a superfluid state as a consequence of the 
spontaneously broken U(l) symmetry. 

There is another way to tackle the slave boson action eq. (|7]). One can decouple the hop- 
ping term by introducing a single Hubbard-Stratonovich field $. The action then becomes 
bilinear in all slave boson fields which can then be integrated out exactly. One then finds 

Z = j D&DX exp Lj P dr(-Y,m]^ + Y,^) ) 

!Jl-exp(-/3(zA 4 + en)) 
where, in the occupation number space, £ n are the eigenvalues of M m>n = E n b m ^ n + 
+ m$5 m+ i jn + y/n$>*5 mjTl+ i, where the E n are the eigenvalues of the atomic problem 
(eq. (0)) and $ is treated in mean- field approximation. Due to the particular form of eq. 
(PI) the constraint can be handled exactly. In terms of the order parameter $, the action 
reads: 

S = -j- |$| 2 -ln(]Texp(-/3a)). (12) 



In the hard core limit, the 2 eigenvalues £ CT are given by 2£ CT = —fi ± y/i 2 + 4 | $ | 2 . This 
result is well known in the context of the ferromagnetic X — Y model and usually serves as 
a pedagogical starting point for the discussion of the paramagnetic-ferromagnetic transition 
in this model fl6| . Solving the saddle-point condition which follows from the action ([12]) 
yields the (/x, T) phase diagram where the phase boundary between the superfluid and the 
normal states is given by \i = 2T c th~ 1 ([i/2dt). This might serve as a description of the 
A-transition in 4 He. Assuming that the A He atoms are sitting on a lattice that is half-filled 
and that they only experience a hard-core interaction yields the ratio T c /Tbe = 0.7, which 
is identical to the ratio T c exp /Tbe as obtained with the help of the experimental data. 



b) FINITE U PROBLEM. 

In the saddle-point approximation the slave-boson action becomes: 

S = ^ I b n | 2 (E n + Ao) — A 



(13) 

n,m 

I minimized numerically both eq. (|1^) and eq. (|13|). Even though both approaches 
yield very different looking equations, it turns out that they deliver identical results at zero 
temperature. The difficult problem of minimizing the slave-boson action is handled in the 
following way. First of all one phase can be removed owing to the gauge symmetry of the 
action. Second the phase of the physical Bose field can not be determined by the saddle- 
points equations, as usual in a superfluid state, but all the others are readily seen to be 
equal to the latter, so as to minimize the kinetic energy. There is thus a single Goldstone 
mode. The constraint expresses the fact that the N bosons to be considered are restricted 
to the surface of a N dimensional hyper-sphere which I parameterize in polar coordinates. I 
am thus left with determining the angles. It turns out that the number of angles that differ 
from zero gradually reduces when the interaction is raised up. This is exemplified in fig. 
1 where I show the amplitude of the first 4 bosonic fields as a function of the interaction 
strength at density p — 1. Even for moderate coupling, say U = t, only the first 5 fields 
differ from 0. This implies that amplitude fluctuations are very substantially weakened as 
compared to the ideal Bose gas. This physical effect persists down to any finite interaction 
strength. 

The Mott insulating state is reached when a single b-field differs from 0. In this case the 
superfluid density, which vanishes, is very different from the density which takes (here) the 
value 1. We meet a very different situation as in the weakly interacting Bose gas theory, 
where both quantities are identical due to Galilean invariance at zero temperature |T^J . How 



they start to deviate from each other is shown in the inset of fig. 1, where I plot | < B > \ 2 /p 
as a function of the interaction strength at density p — 1. 

The phase diagram (fig. 2) shows Mott insulating states corresponding to commensurate 
densities, which are identical to those obtained by ||. They are in good agreement with 
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QMC results by Trivedi and Ullah ||, even though the insulating lobes are somewhat too 
small. This is a drawback of the method which is well-known from similar calculations for 
interacting fermions fli~2" |. This originates in the fact that the action (^|), even though exact, 
does not yield the correct non-interacting limit at mean-field level. Phase fluctuations in 
this approach are somewhat atypical. At finite temperature they are expected to cause 
the superfluid-normal transition. Even though there are a lot of phases that fluctuate, 
only a single one is relevant, all other being massive. This is in contrast to the fermionic 
Hubbard model where all phases but one can be gauged away and the fluctuations of the 
remaining one leads to a massive mode that is then not expected to destroy the condensate 



15| . Here the fluctuations bring a rich excitation spectrum. Considering as an example the 



gaussian fluctuations in the n = 1 insulating lobe leads to a decoupling of the propagator 
matrix. After having integrated out the (real) b\ field and the constraint field and introduced 
ot k = Eq — Ei + |&i| 2 £fc and (3k = E 2 — E\ + 2|5 1 | 2 t fc one obtains: 



uf — ^2 K^i—i^n + E n — E^bn^ 

n>3,k 



( 



\ ( b \ 



(14) 



h* 

y °2,-k j 



-iu n + a k ^/2\bi\H k 

It follows that the spectrum is split into a set of localized high energy excitations (for n > 3) 
and a continuum which arises from the small n part of the action. Looking for a vanishing 
gap provides expressions for the superfluid-insulator lines for the iV-th lobe (N > 1): 



E 

u 



2N - 1 



dt 
U 



± 



(15) 



This is the analytical expression for the Mott-insulating lobes which can be obtained either 
with this method or with the collective field approach. In turn physical response functions 
such as the density-density correlation function can be computed. The Mott gap following 
from the low energy part vanishes at the tip of the lobe where the spectrum is changing 
from being gapful and massful to gapless and massless. At this particular point I calculated 
numerically for N = 1 the density-density correlation function N(q, u) on the 2-d square 



lattice which is displayed in fig. 3. The latter is vanishing identically for k = reflecting the 
incompressibility of the system. As a result charge fluctuations are mostly high energetic 
and inhomogeneous up to a fraction of low but finite energy excitations resulting into a 
gapless incompressible state. The spectrum is very broad and is extending well over the 
band width of the corresponding non-interacting Fermi system. It mostly consists of a two- 
peak structure following from the two modes of the fluctuation matrix. At low momentum 
and energy the imaginary part of N(q, to) is obtained as: 

t ATI \ ©(<*>- eg) q 2 na s 



with the sound velocity c = At/y Q\/2 — 8. This response function thus exhibits an integrable 
singularity at the threshhold. As a result the DC conductivity is finite and takes the universal 
value e 2 /16h where e denotes the charge of the bosons. However a finite value for the 
DC conductivity heavily relies on the particular form of ImN(q,u) as found in eq (|TE). 



Keeping in mind that all self-energy corrections to the slave boson propagators are neglected 
makes it unlikely that such a peaked behavior is a genuine feature of the model. Moreover 
an additional symmetry of the Hamiltonian appearing at a discrete set of points of the 
phase diagram corresponding to the tip of the lobes could not be identified. Thus self- 
energy corrections must exist. Here the holon and doublon propagators only coincide for 
small energy and momentum. The differences are responsible for the 2-peak structure in 
ImN(q,u}) rather than a 1-peak structure. Calculations of N(q,u) away from the tip of the 
lobe yields similar looking results apart from a gap as obtained from eq (|L5|). The detailed 
expression for N(q, u) as well as additional calculations in the superfluid domain that are 



in progress, will be published elsewhere [IS 



Conclusion 

In this paper the superfluid-insulator transition that occurs in interacting bosons systems 
is considered. To this aim an auxiliary boson representation is introduced and I showed that 
the slave boson mean-field theory is identical to a mean-field theory on an order parameter. 
Despite of its apparent complexity it allowed to obtain an analytical expression for the Mott 
insulating lobes as well as for the density- density correlation function. This newly introduced 
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framework is used to show that many energy scales appear in the excitation spectrum. In 
the strong coupling regime it consists of a continuum of low energy excitations and a set of 
localized high energy excitations. The latter are weakening the amplitude fluctuations in a 
substantial way, even for moderate coupling. 
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FIGURE CAPTIONS. FIG. 1 Amplitude of the slave boson fields b n as functions of the 
interaction strength at the commensurate density p — 1. The curves A,B,C,D correspond 
respectively to n — 0, 1,2,3. Inset: | < B > | 2 as a function of the interaction strength at 
the commensurate density p — 1. 

FIG. 2 Chemical potential versus hopping phase diagram at zero temperature. The 
insulating lobes correspond to the densities 1 and 2. 

FIG. 3 Density-density correlation function at the tip of the first Mott lobe on the square 
lattice as functions of frequency. The curves A, B, C, D are calculated for wave-vectors on 
the diagonal of the first Brillouin zone for respectively q x = (1, 2, 3, 
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